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Using a lattice-gas model with pairwise interactions, we study the ordered structures, coverage 
dependence of the heat of adsorption, and other experimentally observable behavior of adsorbed 
CO overlayers on Pd(lOO) single crystal surfaces. Transfer matrix and Monte Carlo methods give 
accurate information regarding the lattice-gas model that often contradicts simple mean-field-like 
analysis. We demonstrate the usefulness of the model by reproducing experimental results over a 
large range of pressures and temperatures. 
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I. INTRODUCTION 



CO adsorption on metal surfaces has been studied 
extensively as a benchmark system for chemisorption. 
In particular, great deal of information has been accu- 
mulated during the last thirty years about CO adsorp- 
tion on Pd(lOO) surfaces using several different experi- 
mental techniques, and under a range of pressures and 
temperatures j ^ i ^ i? There are also detailed theoretical 
studies of CO adsorption on Pd(lOO) using first-principles 
approaches. ^'^ However, surprisingly, no extensive statis- 
tical mechanics studies have been performed to precisely 
determine adlayer ordering and phase transitions for this 
system. This is a significant omission since such analysis 
provides strong constraints on the type and magnitude 
of adspecies interactions. It also provides a reliable de- 
termination of thermodynamic quantities, which is not 
possible with simplified analyses. Such information is in- 
valuable in interpreting other experiments, e.g., related 
to CO adsorption energies. 

A long-standing motivation for such detailed studies of 
simple chemisorption systems is to provide insight into 
catalytic surface reactions. It is also well-recognized that 
ordering and islanding of reactants will limit the utility 
or validity of mean-field type rate equation treatments of 
the reaction kinetics.^ Hence, accurate and robust atom- 
istic modeling of adlayer structure for individual reac- 
tants is a crucial first step in building realistic atomistic 
model of related surface reactions, e.g., for CO oxidation 
on Pd(lOO).^ 

In this paper, we develop and analyze a lattice-gas 
(LG) model for CO adlayers on Pd(lOO). Our focus is on 
equilibrium aspects of this system, since molecularly ad- 
sorbed CO can diffuse quite rapidly on this surface under 
normal situations facilitating adlayer equilibration. How- 
ever, some nonequilibrium issues are also addressed. The 
main techniques used to analyze behavior of the lattice- 
gas model are the transfer matrix method and Monte 
Carlo simulation. We aim to reproduce as many exper- 
imental observations as possible using a relatively sim- 
ple model. We present results regarding surface ordering 
below 0.5 monolayers (ML) in Sec. IIIII the structure of 
dense CO adlayers in Sec. IIVI the heat of adsorption in 
Sec. and the adsorption isobars in Sec. IVII 




FIG. 1: Schematic of the c(2V2 x y/2)M^° ordered structure 
and of the pairwise interactions used in this paper. The large 
open circles represent Pd atoms and the small solid circle 
represents CO molecules adsorbed on the bridge site of the 
Pd(lOO) surface. 



II. LATTICE-GAS MODEL FOR CO/PD(100) 
AND ITS ANALYSIS 



Various experiments^*^ show that adsorbed CO re- 
sides only at bridge sites on Pd(lOO). Thus, our model- 
ing of equilibrated adlayer configurations allows popula- 
tion of bridges sites only. We note, however, that it was 
suggested^ that during adsorption, CO is first steered to- 
wards the less favorable top sites. Thus, in more general 
modeling of nonequilibrium configurations under reac- 
tion conditions, it is appropriate to allow population of 
other sites. ^ Below, a = 2.75 Adenotes the unit cell size 
of Pd(lOO) surface. 

Our LG modeling also assumes only pairwise- 
interactions between CO adsorbates. Figure illustrates 
the specific interactions used: we incorporate nearest- 
neighbor (NN) interactions uji for CO pairs separated by 
distance a/\/2, second NN (2NN) interactions ui2 for sep- 
aration a, third NN (3NN) interactions for separation 
\/2a, and sometimes fourth NN (4NN) interactions W4 
for separation ^hjla. Also illustrated in the figure is 
the experimentally observed c(2V2 x \/2)R\h° ordered 
structure. 

In applying our model to analyze the heat of adsorption 
and the adsorption isobars for CO/Pd(100), we need also 
information of the initial heat of adsorption at low cov- 
erage, which corresponds to the absorption energy of an 
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isolated CO molecule at a bridge site. This quantity was 
measured at 1.55 eV in an early experiment by Tracy and 
Palmbergfi and at 1.67 eV in more recent experimentsi^ 
Using density functional theory (DFT) incorporating the 
generalized gradient approximation (GGA), Eichlcr and 
Hafnes^ studied the potential energy surface for CO ad- 
sorption on Pd(lOO). They reported the adsorption en- 
ergy for CO on a bridge site to 1.92 eV. This value is 
larger than the experimental estimate, but the trend that 
bridge sites are favored over top and hollow sites is con- 
sistent with experimental findings. 

To analyze the above two-dimensional LG model, in 
this paper we use two standard yet powerful statistical 
mechanical techniques: the transfer matrix (TM) and 
the Monte Carlo (MC) methods. These two methods 
are often complementary. Using TM, one can always 
obtain the equilibrium free energy and other thermody- 
namic properties of the system. However, it is difficult to 
include long-range interactions using TM. For example, 
performing analysis on a system or strip of size M x oo, 
in order to include 3NN interaction tJs, it is necessary to 
consider all configurations of two columns spanning the 
strip. Thus, one must consider a total of 2"^^ configura- 
tions, if one does not reduce this number by exclusion and 
symmetry properties. On the other hand, it is relatively 
straight-forward to include long-range interactions using 
the MC method. However, standard MC can become 
inefficient, especially when there are strong repulsive ad- 
species interactions, or for low temperatures. 



III. c{2V2 X y2)i?45° ORDERING BELOW 0.5 ML 

From the observation of c(2\/2 x \/2)R45° ordered 
structure, it has been deduced that lateral interactions 
between CO(ads) consist of strong NN and 2NN repul- 
sions, and a weak 3NN repulsion. It has also been pointed 
out by Behm et al^ that the 3NN repulsion lo^ should 
not be too strong, otherwise it would instead produce a 
(•\/5/2 X ^5/2)i?18.4° structure, which is not observed 
experimentally. 

In order to quantify the effect of the 3NN interaction 
on the ordering of CO adsorbates, we conduct a study of 
the phase transitions of the lattice-gas model with very 
strong NN and 2NN repulsions (i.e., exclusion), and a 
finite 3NN repulsion (wa > 0). Different aspects of this 
model has been studied previously, but the complete pic- 
ture is not available. 

In the case of no 3NN interactions (^3 = 0), it was 
found that there exists a very "weak" transition from a 
disordered phase to a semi-ordered phase upon increas- 
ing the CO coverage above about 6*00 = 0.477.;^- At 
the maximum ^co — 0.5, the semi-ordered phase con- 
sists of alternating half-filled diagonal rows of bridge sites 
which can slide with respect to each other without en- 
ergy penalty. In the limit of very strong 3NN repulsion 
(ws = 00), it was found"'^^ that there is a first-order liquid- 
solid-like transition with increasing CO coverage. The 



ordered phase has a (•\/5/2 x -\/5/2)i?18.4° structure. 

(Traditionally, results of the above studies are de- 
scribed with respect to the square lattice of bridge sites. 
This lattice is rotated by 45° from the square lattice of 
Pd(lOO) substrate atoms, and has twice the number of 
sites as Pd atoms. Thus, the coverage on this lattice sat- 
isfies p = O.50CO- For ^3 = 0, the semi-ordered phase is 
commonly denoted as (2x1) phase, although it has long- 
range order in one dimension only. It has a maximal 
coverage of p = 0.25. For oj^ — 00, the ordered phase is 
described as (a/S x ^/b)R2Q.(!>° with respect to the square 
lattice of bridge sites.) 

Our TM and MC study shows that as uj^ increases 
from zero, the "weak" transition to the semi-ordered 
phase becomes stronger, and for larger 3NN repulsion, 
the transition converts to a first-order transition to 
the c(2\/2 X \/2)-R45° ordered structure. This implies 
the existence of a tricritical point. From MC simula- 
tions, we estimate that this tricritical point is located at 
{(3pl\(3ujI) = (5.2,0.3), where [3 = 1/(A:bT) and n de- 
notes the chemical potential for the adsorbed CO (which 
is described in more detail in the following sections) . The 
corresponding coverage is 9co = 0.46 (or p = 0.23). 

As (3uJz further increases above approximately 1.9, 
there is another first-order transition, but this time, the 
transition is from the disordered phase to a (^^5/2 x 
^5/2)i?18.4° phase. Figure [5] shows our preliminary re- 
sults directed towards mapping out the /i-/3 phase di- 
agram of the model with NN and 2NN exclusion, and 
3NN repulsion. It shows two first-order liquid-solid-like 
transitions (dashed lines) from disordered phase to ei- 
ther a c(2V2 X -\/2)i?45° ordered phase, or a (■y/5/2 x 
^5/2)i?18.4° ordered phase. Also there is a sohd-solid- 
like transition from the (■\/5/2 x y^5/2)i?18.4° phase to 
the c{2^/2 X V2)i?45° phase. There is a tricritical point 
connecting the continuous and first-order transitions to 
the c(2\f2 X \/2)-R45° phase. However, how the termina- 
tion of the (■\/5/2 X ^5/2)_R18.4° transition line is not 
clearly determined from the existing analysis. It is likely 
that the transition line bends towards and merges with 
the c(2\f2 X \/2)^45° transition line sharply at around 
/3W3 = 1.9. 

Under the framework of the lattice gas model with only 
repulsive interactions between neighboring pairs up to 
3NN, one can concl ude tha t loz < 0.5 eV so that at room 
temperature, no (■\/5/2x ^5/2)i?18.4° phase shall be ob- 
served. A caveat is that introducing further neighboring 
interactions can change the phase diagram significantly 
and the above constraint on the magnitude of ^3 is no 
longer valid if longer-ranged interactions should be con- 
sidered. 

Based on the above analysis (and also our investi- 
gations in subsequent sections), we assign a value of 
u!3 = 0.03 eV for the strength of the 3NN repulsive in- 
teractions. In the following analyses which include con- 
sideration of behavior for CO coverage above 0.5 ML, it 
is necessary to relax the constraint of 2NN exclusions. 
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FIG. 2: /i-6' phase diagram of the lattice-gas model with NN 
and 2NN exclusion and 3NN repulsion. The solid line denotes 
a continuous transition obtained from transfer matrix (TM) 
finite-size-scaling^^ (10-12 scahng). Also shown are two first- 
order transitions obtained from locating local maximum in 
dO/d^ using two methods: the dashed lines are from TM (with 
strip of size 10) calculations, and the symbols are from Monte 
Carlo simulations with the histogram method.— 



The value of lj2 of around 0.17 eV is determined from 
our analysis of the heat of adsorption in Sec. and ad- 
sorption isobars in Sec. IVII However, before presenting 
these analyses, in Sec. IIVI we provide a more complete 
picture of adlayer ordering by describing behavior at cov- 
erage above 0.5 ML (using the parameter choice wi = oo, 
U02 = 0.17 eV, and = 0.03 eV.) 
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FIG. 3: Snapshot of Monte Carlo simulations using the 
Glauber dynamics and the Metropolis algorithm for the LG 
model at fixed pressure while lowering the temperature. The 
annealing rate is 10* MCS/K (each site in the system is sam- 
pled once on average for each MCS). Other parameters are 
PCD = lO"'^ Torr, ui ^ 00, L02 = 0.17 eV, ujs = 0.03 eV. 
The snapshot is taken at T = 400 K. For illustration, we de- 
note CO with exactly two 2NN and four 3NN by a circle, and 
all other CO ("defects") by a black dot. Also note that the 
Pd(lOO) substrate (not shown) is rotated 45°. Shown in the 
figure is a L = 128 subsystem in a simulation with L = 256 
using periodic boundary conditions, ^co = 0.511 ML. 



IV. STRUCTURE OF DENSE CO ADLAYERS 

Because of the difficulties of experimental techniques 
(e.g., work function and infrared analysis) in dealing with 
high CO coverage (^co > 0.5), structures of dense CO 
adlayers on Pd(lOO) are a matter of some debate (i^'i^'i^ 
As 9 CO increases above 0.5 ML, it is concluded from 
diffraction studies that the adlayers structure undergoes 
a commensurate-incommensurate transition (CIT). The 
study by Schuster et al^ suggests it is in the Pokrovsky- 
Talapov universality class, consistent with expectation 
from symmetry arguments. 

Lattice-gas models are not ideal for study of CIT's, 
since they put too many constraints on the structure of 
domain walls. Nonetheless, we perform a Monte Carlo 
study to investigate CO adlayer structure above 1/2 ML. 
We simulate the system at a fixed pressure while low- 
ering the temperature. We use Glauber dynamics (cor- 
responding to adsorption/desorption in the LG model) 
with the Metropolis algorithm. Occasionally, we also mix 
in Kawasaki dynamics (corresponding to diffusion in the 
LG model) with the Glauber dynamics. Here, our fo- 
cus is in the equilibrium structure of CO adlayer, thus 
we do not need to mimic the physical kinetics accurately. 



The primary challenge is that due to CO adspecies repul- 
sions, adlayers can become nearly "frozen" using normal 
dynamics at high CO coverage and low substrate temper- 
ature. In fact, the LG model with NN and 2NN exclusion 
has been used to study the glass transitionpii 

Although our simulations are not faithful to the phys- 
ical adlayer dynamics, they provide at least some quali- 
tative insights into adlayer structure. For a fixed system 
size with periodic boundary conditions, if we lower the 
temperature very slowly, eventually a single domain oc- 
cupies the whole system. Upon further lowering the tem- 
perature, defects are formed. However, we are unable to 
observe any well-defined domain wall structure as sug- 
gested by Berndt and Bradshaw.^^ On the other hand, 
if we lower the temperature more quickly, then differ- 
ent domains still occupy the system near the transition 
point. Further decreasing the temperature is accompa- 
nied by the enlargement of those original domains, and 
emergence of defects inside different domains. We show 
in Fig. 13 a snapshot of such a configuration generated by 
Monte Carlo simulations. 
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V. HEAT OF ADSORPTION 

Assuming equilibrium between CO in the gas phase 
and the chemisorbed phase, the Clausius-Clapeyron 
equation relates the gas phase pressure, P, to an isos- 
teric heat of adsorption. Est, via 



[d\nP/d{l/T)]e = -Est/kB. 



(1) 



Various experimental techniquesAiiiSiiS give similar re- 
sults for the value of Est as Oqo ~^ 0, while conflicting 
results have been obtained for the coverage dependence 
of Est- Most studiesiiiSiiS show a decrease in ii^st with 
increasing ^co i while the study by Behm et al^ shows a 
roughly constant Est for 9co < 0.5. Presence of carbon 
is suggested^ as the reason for this discrepancy, a claim 
disputed by others .J^ 

The presence of lateral interactions between CO ad- 
species is often invokedi*i^ to explain the coverage de- 
pendence of -Est. Yet to our knowledge, no systematic 
study of the isosteric heat of adsorption in an interact- 
ing lattice-gas model for this system has been performed 
previously. There are some general studies of effects of 
lateral interactions on the heat of adsomtion using mean- 
field approximations (see, e.g., Ref. |20| and references 
therein) which are not reliable for models with relatively 
short-range interactions. 

In lattice-gas modeling, it is appropriate to use the 
grand canonical ensemble. We assume simply that the 
gas phase pressure P is related to the chemical poten- 
tial through /i = kBT\n{P/ Pq). More accurate forms of 
the relationship between the pressure and the chemical 
potential will introduce corrections to the heat of adsorp- 
tion on the order of /csT, which is negligible for present 
purposes. Unlike the isobar experiments (treated later in 
Sec. IVI|I . for our analysis here, we do not need informa- 
tion regarding prefactors for desorption, or the sticking 
coefficient for adsorption. 

Using the transfer matrix method, we examine the cov- 
erage dependence of Est for the lattice-gas model with 
various choices of interactions. Any coverage dependence 
would reflect the influence of adspecies interactions which 
can cause Est to deviate from its limiting value of eb for 
low coverage. Some of the results are shown in Fig. ^ 

For the lattice gas model with only 2NN repulsive in- 
teraction only, the heat of adsorption is effectively a step 
fimction, with the form 



Es, 



lie < 0.5 



ef, - 4w2 if 0.5 < 61 < 1. 



(2) 



The result can be explained as follows: for < 0.5, ad- 
sorbed CO molecules can rearrange themselves to avoid 
any 2NN pairs. Around 9 — 0.5, they form a near-perfect 
0(2^2 X \/2)i?45° adlayer. To adsorb more CO molecules 
beyond this near-perfect overlayer, it is necessary to cre- 
ate four 2nd NN pairs. See Fig.|Slfor an illustration. 

With longer-ranged repulsive interactions, Est de- 
creases as 9co increases even for 6*00 < 0.5. However, 
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FIG. 4: Transfer matrix calculation of the isosteric heat of 
adsorption versus CO coverage on Pd(lOO). The lines are 
obtained using strips of size 8, and the symbols are obtained 
using strips of size 10. Nearest-neighbor interaction uji is 
assumed to be infinite and other neighboring interactions are 
shown in the figure. T — 300 K. 




FIG. 5: Schematic showing the accommodation of an extra 
CO molecule in an otherwise perfect c(2\/2 x \/2)^45° struc- 
ture. 



the coverage dependence is quite nonlinear. For example, 
with 3NN interactions only (the dotted line) , Est only de- 
creases slightly for 6co < 0.25. This is again due to the 
fact that below this coverage, CO (ads) can easily arrange 
themselves in a way to avoid any 3NN pairs. 

With the present set of parameters for lateral interac- 
tions, the lattice-gas model produces a coverage depen- 
dence of the heat of adsorption somewhat between the 
experimental results of Behm et al. and other groups. 
Most of the decrease in the heat of adsorption occurs 
when 6*00 > 0.5, while only a slight decrease occurs when 
9co < 0.5. 

The near parabolic decrease in Est starting from ^co = 
0, as well as the transient increase after the ordering tran- 
sition, shown by Fig. 0] are quite reminiscent to the ex- 
perimental results of Guo and Yates^i for CO adsorption 
on Pd(lll). However, they reported a plateau at an or- 
dering transition, while we see first a sharp drop and then 
a plateau at the ordering transition point. 
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FIG. 6: Adsorption isobars calculated from the lattice-gas 
model. Two sets of parameters are used, pco ranges from 
10"^ to 0.1 Torr. 



VI. ADSORPTION ISOBARS 

For a CO adlayer in equilibrium with gas phase CO, 
it is clear that as the surface temperature increases (at 
fixed pressure) , the CO coverage will decrease. For higher 
(fixed) pressures, this decrease will be delayed until a 
higher temperature range. To quantify this behavior us- 
ing our LG model, one needs a more quantitative deter- 
mination of the relationship between pressure and chemi- 
cal potential than that presented in the previous section. 
To this end additional assumptions are needed. We as- 
sume that the impingement rate is given by P/v^TrmfegT 
and the attempt frequency for desorption is i^q, then we 
assume that 



Note that we "naively" choose vq 



10 



13 , 



for the 



h'Q\/2TrmkBT 



(3) 



The initial sticking coefficient when ^co = is taken to 
be unity. Thus by assuming equilibrium of CO between 
gas phase and the chemisorption phase, one can calculate 
the adsorption isobar of the lattice gas model using either 
the transfer matrix or the Monte Carlo method. 

The results are shown in Fig. |H| For low pressures, as 
T decreases, 9co first increases, then reach a plateau at 
^co = 0-5- For Pco — 10^^ Torr, the plateau occurs 
near T = 420 K. This result is in very good agreement 
with experimentsi^ii^ Upon a further decrease in temper- 
ature, ^co again increases above 0.5 ML. Experimentally, 
this corresponds to the adlayer moving into the regime of 
the commensurate-incommensurate transition discussed 
in Sec. IIVI Here, we simply note that the temperature 
where this occurs depends sensitively on the value of the 
2NN interactions uj2- By choosing LU2 = 0.17 eV, the tran- 
sition occurs between 340 K to 400 K for pco between 
10-9 to 10"'^ Torr, in agreement with experiments^Siia 
It is also significant that at high pressures, the LG model 
predicts disappearance of the plateau near 6*00 = 0.5, 
which is also observed experimentallypi^ 



desorption prefactor, while Behm et al. obtain a value on 
the order of 10^^ s^^ from their adsorption isobars. Con- 
sequently there is some discrepancy between the lattice- 
gas model prediction and experiments at low coverage. 
Specifically for 6co < 0.5 ML, significant desorption oc- 
curs at a higher temperature than experiments. Adopt- 
ing adsorption energy e;, from DFT calculations and as- 
suming the same prefactor makes the discrepancy even 
larger. Also the discrepancy will not likely be resolved by 
modification of interactions, since at low coverage lateral 
interactions are quite insignificant. 



VII. SUMMARY 

We have performed a combined transfer matrix and 
Monte Carlo study of a lattice-gas model for CO adlay- 
ers on Pd(lOO). Model predictions are compared against 
a variety of experimental observations. Of particular sig- 
nificance is our determination of repulsive adspecies in- 
teractions: LUi = 00, W2 = 0.17 eV, and W3 = 0.03 eV. 
Our estimate of the 2NN interaction of 0.17 eV agrees 
well with Wu and Metiu.^ Our assignment of a weak 3NN 
interaction is consistent with early qualitative arguments 
by Behm et al.^ Our use of a binding energy of eh — 1.6 
eV in analysis of the heat of adsorption and the adsorp- 
tion isobars is consistent with experimental estimates, 
but this value is somewhat smaller than that from DFT 
which is close to 2 eV. 
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